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One of the most distinctive features of graphene is its huge inter-Landau-level splitting in experi¬ 
mentally attainable magnetic fields which results in the room-temperature quantum Hall effect. In 
this paper we calculate the longitudinal conductivity induced by two-phonon scattering in graphene 
in a quantizing magnetic field at elevated temperatures. It is concluded that the purely phonon- 
induced scattering, negligible for conventional semiconductor heterostructures under quantum Hall 
conditions, becomes comparable to the disorder-induced contribution to the dissipative conductivity 
of graphene in the quantum Hall regime. 

PACS numbers: 73.43.Cd, 73.22.Pr 


I. INTRODUCTION 

The quantum Hall effect, discovered in 1980 by Klaus von Klitzing^. allows one to determine the quantum resistance 
standard in terms of the electron charge and Plank’s constant with a spectacular accuracy. However, the level of 
precision necessary for metrology applications (a few parts per billion) requires for conventional quasi-two-dimensional 
semiconductor systems the use of ultra-low temperatures and high magnetic fields/S.^ The discovery of graphene^ in 
2004 led to a revival of interest in quantum Hall effect physics—£ The energy of Landau Levels (LLs) in graphene is 
given by 

En = ±(Hvf/Ib)V2N, ( 1 ) 

where up = 1 x 10 6 m/s is the Fermi velocity, Ib = y/h/eB is the magnetic length, B is the magnitude of the 
applied magnetic field, N is the LL number, and the “±” signs correspond to the electrons and holes respectively— 
For B = 10 T the separation between the zero and the first LL in graphene is AKoi/Zcb ~ 1300 K (in comparison, 
AE/Rb ~ 200 K in conventional quasi-two-dimensional semiconductor systems). Therefore, the quantum Hall effect 
in graphene can be observed even at room temperatures in experimentally attainable magnetic fields. Indeed, the 
room-temperature quantum Hall effect in graphene was first observed in 2007^2 however, extremely high magnetic 
fields (up to 45 T) were needed. 

It has been understood since the 1930’s that the longitudinal conductivity of metals in a quantizing magnetic field 
increases with increasing electron scattering:— To develop a graphene-based quantum Hall standard of resistance^ — 
which would work at elevated temperatures and in moderate magnetic fields, it is essential to examine the contri¬ 
butions from different scattering processes to the longitudinal conductivity a xx in the quantum Hall regime (QHR), 
T < AE m /kB, as a xx provides the major correction to the quantized value of the Hall resistance. In conventional 
semiconductor systems the quantum Hall effect is observed at liquid helium temperatures only, and the value of a xx 
is governed by scattering on disorder. In the QHR, <j xx depends exponentially on the energy separation between 
the Fermi level and the nearest LL. The longitudinal conductivity prefactor has a theoretically-predicted value of 
e 2 //i in the presence of short-range disorder and 2 e 2 /h in the presence of long-range disorder.—^ It was shown that 
the contribution from phonon scattering is negligibly small. The exceptions are exotic cases such as magnetoroton 
dissociation, in which the whole effect arises from scattering on phonons^ or an enhancement of phonon-induced 
scattering near the intersection of two LLs corresponding to different size-quantization subbands—*— 1 However, as 
we show in this paper, the plronon-scattering mechanism dominates in the high-temperature QHR in graphene, since 
at T > Ti b = [fhs/l b)/ kB (where s = 2 x 10 4 m/s is the sound velocity in graphene 22 ’ 23 ) the energy of an acoustic 
phonon with a wavevector comparable to the inverse magnetic length is much smaller than the temperature; therefore, 
the number of such phonons increases drastically. 

In this paper we restrict our consideration to electron scattering induced by the interaction with intrinsic in-plane 
phonons only, neglecting the effects of disorder as well as electron-electron interactions. Our calculations result in the 
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lower estimate of the longitudinal conductivity a xx in graphene in the QHR, as we do not study electron interactions 
with various other types of phonons such as out-of-plane flexural phonons, which are present in suspended samples; 24 ’ 25 
or bulk acoustic phonons in a substrate, which interact with electrons in graphene on polar substrates such as boron 
nitride^ or silicon carbide via piezoelectric coupling —L2& 


II. THEORETICAL FORMALISM: TWO-PHONON SCATTERING WITHIN ONE LANDAU LEVEL 


The wave function of a charged carrier (electron or hole) in graphene subjected to a magnetic field B normal to 
the graphene plane is given in the Landau gauge, A = (0,I3x,0), by the following expression^— 




C N 


N,ky — 




exp 


1 

IT 

i 

h 

o 

to 

_1 

H n [(a: - x 0 ) /l b\ 

1 

CN 

1 

_l 

±V2NH N - 1 [(x — xq ) / Ib ]. 


exp ( ik v y ) 

\[Ly 


( 2 ) 


where Co = 1 and Cn^o = l/"\/2, Xo = l^k y is the guiding centre coordinate, k y is the electron wavevector y- 
component, H jy are the Hermite polynomials, L x and L y are the graphene sample dimensions, the “±” sign corre¬ 
sponds to electrons and holes, and the spin and valley indeces are omitted. Here the z-axis has been chosen in the 
direction of the applied magnetic field and the x-axis and y-axis are in the graphene plane. The wave function given 
by Eq. @ is defined for one Dirac point only and its two components correspond to graphene’s two sublattices. In 
what follows all results will be obtained for electrons only. Calculations for the non-equivalent Dirac point and holes 
can be easily repeated in a similar fashion. From here on we will omit the “±” sign in order to simplify notation. We 
account for the valley and spin degeneracy by multiplying the final result for a xx by the factor of four at the final 
stage of calculations. 

Electron scattering on phonons leads to a change in the k y component of the electron wavevector, which results in 
a change to the electron guiding centre coordinate xq. Electron transitions between LLs due to one phonon scattering 
are suppressed due to the large energy gaps between LLs in graphene in the QHR. Indeed, the number of phonons 
with the energy required for such transitions is n q ~ exp (— AE^^/kBT) (AE^ 1 n 2 is the energy gap between two 
different LLs in graphene), which is very small in the QHR. Inter-LL scattering on acoustic phonons has an additional 
exponentially strong suppression in the matrix element of transition which is markedly different from the case of 
conventional semiconductor systems as discussed in the Appendix. It is also evident that in the QHR only scattering 
on acoustic phonons can provide a noticeable contribution to the longitudinal conductivity <j xx . The optical phonon 
energy in graphene corresponds to the temperature range 1800 K — 2300 K 29 i 30 which leads to very small optical 
phonon occupation numbers at room temperatures and below. In this work, we are interested in phonon-induced 
equilibrium longitudinal conductivity, whereas, magneto-phonon resonance associated with optical phonons has been 
studied in Ref— ir 33 . One-phonon scattering within the same LL is ineffective because of the small width of the LLs, 
which can be arguably achieved in pristine graphene samples. Interestingly, even in heavily disordered samples the 
broadening of LLs was found to shrink in several particular cases For our calculations we assume an infinitely 
narrow band of delocalized states in the middle of each LL with a vanishing density of states in between LLs. In 
this limit one-phonon scattering within one LL is forbidden by energy and momentum conservation. The next order 
process to consider is two-phonon scattering through virtual states with no change in the electron initial and final 
LL numbers but with the change in its in-plane momentum (or the guiding centre coordinate). We calculate the 
longitudinal conductivity cr xx at the Nth LL due to two-phonon scattering using the generalized form of the Einstein 
relationi 21 i 36 i 37 
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where v N = {exp[(£)v — E F ) /k B T] + l}^ 1 is the LL filling factor, E N is the LL energy defined by Eq. ([T]), E F is 
the Fermi energy, k y and k' y are the y-components of the electron wavevector before and after scattering, D is the 
diffusion coefficient, and Wk v ^>k' y is the probability of scattering. 

Two-phonon scattering in graphene in the QHR is possible through two different virtual intermediate states: a 
phonon with wavevector q + is first emitted or a phonon with wavevector q _ is first absorbed. Transitions changing 
the electron LL number in the intermediate states are strongly suppressed due to small values of corresponding matrix 
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elements and the presence of large denominators in the expression for W ky ^, k ' (see Appendix). Therefore, in what 
follows we consider only transitions with no change in the electron LL number. Then, the probability of two-phonon 
scattering is given by Fermi’s golden rule with the matrix element 


M k y ,k’ y (q >q + ) = n q - ( n q + + l) x ^2 
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where n q = [exp (hflg/ksT) — 1 ] 1 is the phonon occupation number, £l q = sq is the acoustic phonon frequency, and 
Vq is the electron-phonon coupling operator. 

Operators describing electron scattering on intrinsic longitudinal (LA) and transverse (TA) acoustic phonons in 
graphene ar o 38 ’ 39 


V LA 

q 


U q q exp (iqr) 


( i( Jd 

\-g h e- i2Lp 


ghe i2ip \ 
igd ) ’ 


( 6 ) 


Vq TA = U q qexp (iqr) ( ig J-i 2 v W q ) . (?) 

where U q = {L x L y )~ lG (h/2pVL q ) 1 ^ 2 , p is the graphene 2D mass density, r is the position vector in the graphene plane, 
and ip is the angle between the phonon wavevector q and the axaxis. The diagonal matrix elements in Eqs. ©-© 
describe electron coupling to the phonon-created deformation potential, and the off-diagonal matrix elements originate 
from the phonon-induced bond-length modulations, which effect hopping amplitudes between two neighbouring sites. 
The corresponding coupling constants were estimated as gd ~ 20 — 30 eV and g k ~ 1.5 — 3.0 eVi 22 i 23,38 ~ — Note 
that the deformation potential couples electrons with LA phonons only. Furthermore, since the second component of 
the electron wave function defined by Eq. © vanishes for the zero LL, electrons in this LL do not interact with TA 
phonons. 

Substituting the electron-phonon scattering operators given by Eqs. ©-© into Eq. © yields 
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where p = {LA,TA}, 7 = {d, h} are the indices introduced to separate the contributions to the longitudinal conduc¬ 
tivity Oxx from LA and TA phonons and from the two different scattering mechanisms discussed above, G is the 
generalized electron-phonon coupling constant with the values Gla,cL = gd , Gla,1i = < 7 /i, GrA,d = 0 , GrA,h = gh , and 
the matrix elements ’ 7 fc , (q) are given by 


M k y A kl (q) = * J *N,k' y ( r ) ex P (±*q r ) ™N,k, (r) dr = iC 2 N exp (±/3 - a) [L% (a) + L° N _ ± (a)] S k ^ ky±qy , (9) 
M k A kl (q) = J ^N,k' v ( r ) ex P (±*q r ) ^LA^N,k y (r) dr 

= - C%N~ 1/2 l B \q x cos 2(p± [k y - k' y ) sin2<^] exp (±/3 - a) L 1 N _ 1 (a) 6 k > k y ±q v , (10) 


M k y A kl (q) = J ®N,k' v ( r ) ex P (±*q r ) ^TA^N,k y (r) dr 

= ~C 2 N N~ l/2 l B [q x sin 2ip =F (k y - k' y ) cos 2 v?] exp (±/3 — a) L 1 N _ 1 {a)5 k ^ ky ± qy . (11) 


In Eqs. ©-(HU), a = (l 2 B / 4) (k' y — k y ) + q 2 , /? = i (l%q x / 2) ( k' y + k y ), N ^ 0, the “±” sign refers to emitted and 

absorbed phonons, L° N and L ] v are the Laguerre polynomials, I is the 2x2 identity matrix, and &la, &ta are given 
by 
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Eq. © is also valid for the zero LL when substituting L 0 N _ X = 0 and L 1 N _ 1 = 0. Note the N 1 ^ 2 factor in Eqs. (ITO - 
m which suppresses the off-diagonal contribution to cr xx in higher LLs. Summation over all possible values of k” in 
Eq. © results in the following expressions for two-phonon scattering matrix elements 
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= ( g h U q ) A [q/ (hs)} 2 N 2 [L 1 n _ 1 (l 2 B q 2 / 2)] 4 exp (- l B q 2 ) 

x 1% (<lZ sin2^ + + q+ cos2 <p + ) 2 (q~ sin 2^" + q~ cos 2(p ~) 2 

x sin 2 [1% (q+ q- - q~q+) /2] S k , ky+q - _ q +, (14) 


MIX C » (9HU g )* [q/ (hs)} 2 N- 2 [L 1 n _ 1 (l 2 B q 2 / 2)} 4 exp (-l 2 B q 2 ) 

X 1% {qt cos2iyj + - q+ Sin2(^ + ) (q~ cos 2 <p~ - q~ sin2<*r) 

x sin 2 [1% (q+q~ - q~q+) /2] 5*, |fcy+ ,-_ 9 +. (15) 

Substituting the calculated probability of the two-phonon scattering 114,, into Eq. © and performing summation 
over k' y as well as integration over all possible values of q + and q yield the following result for the longitudinal 
conductivity 

_ _ (~LA,d I ~LA,h ~TA,h\ /i \ 

Vxx — \& xx + a xx + a xx JVn^—Vn), 
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Here Jo is the Bessel function of the first kind. The expressions for ai x '~ f were multiplied by a factor of four to account 
for the valley and spin degeneracy. From Eqs. G9-GB it is evident that a xx ’ d oZ. x ^ TA ' A due to the relatively 
small value of gh comparing to ff ^ 22]23 ' 38 ~ — and the small N~ 2 factor contained in Eq. m- 

There are two distinctive temperature limits. In the low-temperature limit, T -C T; s , the main value of the integrals 
in Eqs. asD-aiD is formed by q < k B T/hs and ql B can be considered as a small parameter. Expanding the integrand 
into power series in ql B and taking into account only the lowest power term yields the following expression for 
in the low-temperature limit 

- (e 2 /h) (2V/15) [g 4 J ( l 2 B p 2 h 2 s 6 )] (T/T lB f , (18) 

tZx /TA ' h - (e 2 /h) [A^/N 2 ] (55287t 11 /1365) [g 4 / ( l 2 B p 2 h 2 s 6 )] (T/T 1b ) 12 . (19) 

. AN—A 

In Eq. (fTSl) the coefficients are defined by \L 1 n _ 1 ( x / 2 )} = YO A^~ 4 x^, where N > 1. For the zeroth LL, 

3 =0 

oZx^ TA ' h = 0. Note that 5 xx ' d in Eq. (1151) does not depend on the LL number. The temperature dependencies 
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given by Eqs. (1181191) are different from the case of conventional semiconductor heterostructures for both two-phonon 
scatterin g 21 ! 36 ! 37 and phonon-assisted hopping conductivity ! 16 ! 17 In Eq. (1T51) . which corresponds to the deformation 
potential scattering mechanism, the lower power in the temperature dependence of mobility compared to that obtained 
in Refs. }2lll36lf37j is due to phonons in graphene being two-dimensional. As was mentioned above, the contribution 
to longitudinal conductivity at higher LLs given by Eq. (fl9l) is a distinctive feature of graphene without analogy in 
semiconductor systems. 

In the more interesting high-temperature limit, T > Ti b , the main value of the integrals in Eqs. m-m is formed 
by q < 1/Ib, and Hsq/kBT can be considered as a small parameter. By expanding the integrand into power series 
in hsq/ksT and taking into account only the lowest power term we obtain the following result for o 7 ^ 7 in the high- 
temperature limit 


« L J' d = {e 2 /h) (C%$ d N /2\) (. g 4 d /l 2 B p 2 h 2 s 6 ) (T/T lB ), 


where 


AN 

= ■ 2 
1=0 




and 


where 


5% ,TA ' h = (e‘/h) N-> [<!*_,/ (2 6 jt)] [g'J (l^ftV)] (T/T ,.), 


AN-A 

*N-1= E <~ lr U + 
1=0 


a - 2 UI 


( 20 ) 

( 21 ) 

( 22 ) 

(23) 


In Eqs. (I2ll) - (l23l) . T is the Gamma function and 2 -P 1 is the hypergeometric function. In Eq. (j22]l . N 7 ^ 0, since in the 
zeroth LL a x A ^ TA ' h = 0. The coefficients A™ 1 are the same as in Eq. (flTfl) and the coefficients Bj 1 are defined by 

, AN 

[L°n ix/2) + L° N _ X (x/2)] 4 = £ Bfx*. 

1=0 


III. DISCUSSION AND CONCLUSION 

In Fig. |T]we plot a xx = ^ x A,d + ^ xx ’ h +^ x A,h obtained by numerical integration of Eqs. (fl6l) - (fT71) from 0 to 300K for 
B = 10 T (T[ b « 22 K). It can be seen from Fig. Oja) that for T , a xx linearly increases with temperature in line 




FIG. 1: (Colour online) (a) The longitudinal conductivity prefactor a xx as a function of temperature for B = 10 T. (b) 
Expanded view in the 0 — 30 K temperature range. 

with Eqs. (EU1) - (B^1) . From Fig.QJb) one can see the change of the a xx temperature dependence from the high-power 
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to the linear law occurring around T = Ti b . Substituting the numerical values of all the constants into Eqs. (l20l) and 
01 results in the following simplified expression for a xx at elevated temperatures 


a xx » a N (T /300 K) (B/10 T) 1/2 , (24) 

where cfjv has the following values for the six lowest LLs: do = 0.65 e 2 /h, d\ = 0.06e 2 /h, &2 = 0.20e 2 /ft., <73 = 0.19 e 2 /h, 
<74 = 0.15e 2 /h, <75 = 0.14 e 2 /h. Here we used the unscreened value of the deformation potential, which is arguably 
appropriate in the strong magnetic field as the electron motion is quantized and the screening is suppressed. 

There is a natural question of the applicability of the lowest allowed (in our case second) order of perturbation 
theory in electron-phonon interaction at elevated temperatures. This problem was studied in detail in Ref. [.37t | for 
conventional semiconductor quantum wells and bulk acoustic phonon scattering. To compare our results with the 
analysis provided in Ref. [37 ] it is necessary to express them in terms of the diffusion coefficient 
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In Ref. [13 the diffusion coefficient for T > T) B in the lowest (second) order of perturbation theory is written as 
D = ( aT/Ti B ) sIb • Thus, there is a simple connection between the dimensionless electron-phonon interaction 
constant a and our dimensionless constants d/v/ (e 2 /h). Namely, 


a = 0.27 




According to Ref. [37] higher orders of perturbation theory can be neglected for T < T c = Ti B /a. Unlike the case 
of conventional semiconductors, for which a oc lg 2 , for graphene a oc Ig 1 ; therefore, the temperature T c defining 
the validity of perturbation theory becomes independent of magnetic field. For the zero Landau level and for the 
parameters used in our calculations, T c ss 220 K. Above this temperature the phonon-induced mobility is expected 
to change from linear to sub-linear temperature dependence and eventually to saturate. For higher Landau levels 
the perturbation theory cut-off temperature T c well exceeds 300 K and the lowest order perturbation analysis is fully 
valid for ambient conditions. 

In conclusion, we obtained the value of the longitudinal conductivity in graphene in the quantum Hall regime 
due to two-phonon scattering at elevated temperatures which is comparable to the disorder-induced longitudinal 
magneto-conductivity in conventional semiconductor heterostructures I 16,17,13 The predicted distinctive temperature 
and magnetic field dependence of the phonon scattering contribution to the pre-exponential factor in a xx given by 
Eq. (1241) can be easily separated from the temperature- and field-independent contribution caused by disorder when 
analysing experimental data. This should allow the parameters of electron-phonon interaction in graphene to be 
extracted with enhanced accuracy. 
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Appendix: Electron-phonon scattering involving different Landau levels 

The probability of two-phonon scattering through a virtual intermediate state is calculated using Fermi’s golden 
rule as given below 
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Here £), E mi and Ef are the energies of the electron LLs in the initial, intermediate and final states. Clearly, the 
transitions changing the electron LL number in the intermediate states are suppressed compared to the transitions 
conserving the LL number due to the presence of large denominators in Eq. (1A.1I) . 

Notably, single-phonon scattering on acoustic phonons with a change in the LL number is also suppressed because 
of the very small value of the corresponding matrix element 
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Due to energy conservation, q = ^ UF (\JNf — y/Ni ), which results in the following estimate 
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Interestingly, unlike the case of a conventional quasi-two-dimensional semiconductor system, 44,45 the exponential factor 
in Eq. (IA.2I) does not depend on the magnetic field. Instead, it has a non-trivial dependence on the difference between 
the LL numbers of the two involved levels and becomes non-vanishing for very high adjacent LLs. However, these high 
levels are not relevant for the high-temperature quantum Hall effect, which is the subject of our interest; whereas, 


for Ni = 0 and Nf = 1, (Nf 
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